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ABSTRACT 

Observations of the Galactic Centre show evidence of one or two disc-like structures 
of very young stars orbiting the central super-massive black hole within a distance of 
a few 0.1 pc. A number of analyses have been carried out to investigate the dynamical 
behaviour and consequences of these discs, including disc thickness and eccentricity 
growth as well as mutual interaction and warping. However, most of these studies 
have neglected the influence of the stellar cusp surrounding the black hole, which is 
believed to be 1 — 2 orders of magnitude more massive than the disc(s). 

By means of iV-body integrations using our bhint code, we study the impact of 
stellar cusps of different compositions. We find that although the presence of a cusp 
does have an important effect on the evolution of an otherwise isolated flat disc, its 
influence on the evolution of disc thickness and warping is rather mild in a two-disc 
configuration. However, we show that the creation of highly eccentric orbits strongly 
depen ds on the gra ininess of the cusp (i.e. the mean and maximum stellar masses): 
While Chang (2009) recently found that full cycles of Kozai resonance are prevented by 
the presence of an analytic cusp, we show that relaxation processes play an important 
role in such highly dense regions and support short-term resonances. We thus find that 
young disc stars on initially circular orbits can achieve high eccentricities by resonant 
effects also in the presence of a cusp of stellar remnants, yielding a mechanism to 
create S-stars and hyper-velocity stars. 

Furthermore, we discuss the underlying initial mass function (IMF) of the young 
stellar discs and find no definite evidence for a non-canonical IMF. 

Key words: black hole physics - stellar dynamics - methods: TV-body simulations - 
Galaxy: centre. 



1 INTRODUCTION 

Observations of the Galactic Centre revealed one or two 
discs of ~ 6Myr old stars orbiting the central super- 
massive black hole (SMBH) at a distance o f ~0.1 pc 
ilLevin fc Beloborodovll2003l; iPaumard et all [20061 : iLu et ail 
l200fj , 120091 : iBartko et al.ll2009l ). Until now, a large number 
of studies have been carried out to investigate the dynamical 
beh aviour of such s tellar discs around SMBHs. 

lOuadra et al.l (|2008l ) have shown that a single cold 
disc in the absence of a perturbing potential cannot ex- 
plain the obser v ed la rge inclinations and eccentricities. 
iNavakshin et al.1 l|2006h have provided upper mass limits 
to the discs' masses by studyin g the amount of warpin g of 
two discs at large inclination. In lLockmann et al.l (|2008l ) we 
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have shown that a significant number of disc stars eventually 
achieve highly eccentric orbits due to Kozai interaction, thus 
providing a mechanism to create the observed S-stars and 
hyper-velocity stars by disruption of close-passage binaries. 

All these simulations have been carried out neglecting 
the stellar cusp pre s ent in t h e vici nity of the SMBH (e.g., 
ISchodel et"aill2007l '). IChand i|2009f ) argued that a smooth 
background potential hinders the creation of high eccen- 
tricities in full Kozai cycl es since it causes orbital preces - 
sion. On the other hand, lLockmann fc Baumgardtl (|2009h 
have shown that relaxation processes in a grainy stellar 
cusp surrounding the SMBH act on time sca les compara- 
ble w ith the age of the stellar discs (see also iPerets et al.l 
2008). They demonstrate that due to warping and relax- 
ation, two highly inclined stellar discs cannot be recognised 
as fiat circular discs after 5 Myr of interaction. Significant 
contributions arise from both two-body relaxation due to 
stellar encounters customary in star clusters, and resonant 
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relaxation effects, an enhanced rate of r elaxation of angular 
mom e ntum in near-Keplerian systems dRauch fc Tremaind 
1 19961 : iHopman fc Alexander! 120061 : lEilon et al.l l2009h . Fur- 
thermore, both disc and cusp stars contribute significantly 
to the disc evolution. 

Considering the apparent importance of the influence 
of a stellar cusp on the disc evolution, in this paper we chal- 
lenge the results of earlier studies by including both ana- 
lytic and iV-body cusp potentials in corresponding models. 
In Section [2] we discuss the structure of the stellar cusp as 
predicted by theory and derived from observations. Section[3] 
describes the numerical method underlying our integrations 
and our models of the stellar discs. We present our results 
with respect to the studies mentioned above in Section [4] 
and discuss them in Section [S] 



2 STRUCTURE OF THE STELLAR CUSP IN 
THE GALACTIC CENTRE 

The discs of young stars orbiting the SMBH in the centre 
of the Milky Way are embedded in a population of old stars 
and stellar remnants. Theoretical arguments and TV-body 
simulations show that a stellar system around an SMBH 
evolves in to a cusp with a 7 = 1.75 power-law density d is- 
tribution ([BahcaU fc Won] 1 19761; fBaumgardt et al.ll2004aH bl; 
iPreto et al.ll2004l ; lAmaro-Seoane et al.ll2004h . The density 
profile of the stellar cusp in the Galactic Centre is esti- 
m ated from the observ ed luminosity profile and kinematics 
bv lSchodel ei~al1 l|2007h as 

Pm (r) = (2.8±1.3)xlO 6 M pc- 3 (^— ) , (1) 

where r is the distance from the SMBH, 7 = 1.2 inside 
0.22 pc and 7 — 1.75 outside 0.22 pcQ Hence, the stellar cusp 
has a mass within 1.5 pc comparable to the black hole mass. 
We use this result (omitting the uncertainty specified) when 
modelling the cusp as an analytic potential. However, as we 
will discuss with the results in Section [3] relaxation between 
stars can play an important role in the highly dense central 
parsec of our Galaxy. Hence it is not sufficient to include 
such a smooth potential: To perform realistic simulations of 
disc dynamics in the Galactic Centre, it is required to con- 
sider stars and stellar remnants in the iV-body calculations. 

A model of the stel lar populat i on in the central 100 pc of 
our Galaxy is given by I Alexander! l|2005l ). Due to dynamical 
friction causing massive stars to sink towards the centre, 
it can be assumed that the central few 0.1 pc, where we 
also find the young stellar discs, are mass-dominated by the 
most massive components, namely stellar-mass black holes 
i|Freitag et alj|2006t [Alexander fc Hopmanll2009l ). 

In this paper we will analyse the dynamics of the stel- 
lar discs, which will be described in more detail in the next 
section. This restricts our interest to the po tential within 
1 pc fr om the SMBH. At 2pc from the centre [Glisten et al.l 
(|l987h found a rotating ring of gas of 10 6 Mr . the so-called 
circum- nuclear disc (CND). ISubr et ah! |200ct ) have shown 
that the torque exerted by the CND may cause (additional) 



warping of a stellar disc at ~ 0.1 pc from the SMBH. Fur- 
thermore, it may cause somewhat higher eccentricities due 
to Kozai interaction (cf. Section I4.4|l . Throughout this pa- 
per, we will ignore the effect of the CND. This leaves the 
stellar cusp as the only dynamical perturbation to the disc 
evolution. An exhaustive study of the dynamics of stars in 
the Galactic Centre will have to include the effects of the 
CND. However, the general conclusions about the influence 
of a stellar cusp on the disc dynamics as discussed in this 
Paper will remain unaffected. 

3 NUMERICAL METHOD AND MODELS 

To test the mutual interaction of two stellar discs orbit- 
ing the central SMBH, we perfor med direct Af-body inte- 
grati ons using our bhint codeQ I Lockmann fc Baumgardt] 
2008). bhint has been developed specifically to calculate the 
dynamics of stars orbiting around an SMBH, avoiding the 
large sec ular error other A^-body c odes based on the Hermite 
scheme (|Makino fc Aarsetblll992T ) accumulate in case of the 
highly eccentric orbits and extreme mass ratios present in 
the Galactic Centre, bhint integrates the equation of mo- 
tion of a star orbiting an SMBH by moving it along the Ke- 
plerian orbit and adding the forces from other components 
(disc stars and cusp black holes) as perturbations. It solves 
Kepler's equation analytically (and thus without secular er- 
ror) to determine the orbital motion around the dominating 
SMBH, which is assumed to rest at the origin of the coor- 
dinate system, while using the Hermite scheme to integrate 
the perturbing forces. 

bhint includes special treatment of close encounters 
and fast approaches to allow for high precision. In addition, 
it comprises post-Newtonian (PN) approximation of the ef- 
fects o f general relativity up to order 2.5 following iBlanchet] 
(2006), which is added to the perturbations. However, since 
the PN terms need relatively high accuracy and thus a sig- 
nificant amount of computational time, we usually switch on 
these terms only in integrations of highly eccentric orbits, 
where pe ricentre-shift becomes important (see Table [TJ. For 
example, lLockmann fc Baumgardt! (|2009l ) have shown that 
in the regime of the stellar discs, pericentre shift due to the 
spherical cusp is at least one order of magnitude more effec- 
tive than the relativistic one, as long as the eccentricity is 
below e w 0.9. 

Unless stated otherwise, our model parameters are as 
follows: We start with two cold stellar discs around a 3.5 x 
10 6 Mq SMBH at high relative inclination i. The discs have 
finite radial extents ri,2, a surface density E(r) cx r~ 2 ' 5 , and 
a thickness (standard deviation) of 1.44°. The stellar orbits 
have a mean eccentricity of 0.03. We vary the initial values 
for disc masses, inclination, and radial extent within ranges 
compatible with the observations (see Table [TJ. 

The underlying IMF in most of our models has a slope 
of a — 1.35, which was derived by IPaumard et all (|2006T l 
for the stellar discs in the Galactic Centre. This is consider- 
ably flatter t han the c anonic al initial mass function (IMF) 
according to iKroupal (|200lh . f (m) oc m ai , with ao = 



1 These results are con sistent with those of a more recent analysis 
of lSchodel etafl <200Sh . 
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Model Mi n M 2 T2 i m a / bin T end M SB H r SBH i PN 



la 


15 000 


0.05-0.5 


- 


- 


- 


1-120 


1.35 


30% 


5Myr 


15 


0.22 


- - 


lb 


15 000 


0.05-0.5 


- 


- 


_ 


1-120 


1.35 


30% 


5Myr 


- 


- 


yes 


lc 


15 000 


0.05-0.5 


- 


- 


- 


1-120 


1.35 


30% 


5Myr 


- 


- 




2a 


10 000 


0.05-0.5 


5 000 


0.07-0.5 


88° 


1-120 


1.35 


30% 


5Myr 


15 


0.22 




ZD 


10 000 


0.05-0.5 


o uuu 


u.u i — u. 


88° 


1-120 


l.OO 


30% 


5Myr 






yes 


2c 


10 000 


0.05-0.5 


5 000 


0.07-0.5 


88° 


1-120 


1.35 


30% 


5Myr 








3a 


18 000 


0.05-0.5 


9 000 


0.07-0.5 


88° 


1-120 


2.3 


40% 


5Myr 


15 


0.22 




3b 


13 000 


0.05-0.5 


8 500 


0.07-0.5 


105° 


1-120 


2.3 


13% 


5Myr 


15 


0.22 




3c 


11000 


0.05-0.5 


6 500 


0.07-0.5 


88° 


1-120 


2.3 


13% 


4Myr 


15 


0.22 




3d 


24 000 


0.05-0.5 


12 000 


0.07-0.5 


88° 


0.01-120 


canonical 


30% 


4Myr 








4a 


20 000 


0.05-0.5 


10 000 


0.07-0.5 


125° 


1-120 


1.35 


30% 


6Myr 


15 


0.22 


if e > 0.9 


4b 


15 000 


0.05-0.5 


7500 


0.07-0.5 


88° 


1-120 


1.35 


40% 


6Myr 


15 


0.22 




4c 


20 000 


0.03-0.4 


8 000 


0.05-0.4 


100° 


1-120 


2.3 


40% 


6Myr 


15 


0.22 




4d 


15 000 


0.05-0.5 


7500 


0.07-0.5 


88° 


1-120 


1.35 


40% 


6Myr 


10 


0.22 


yes 


Table 1. 


Initial parameters 


used in our 


models. 


These columns give 


the following 


in order: Model 


name; 


mass of the first disc; radial 



extent of the first disc; mass of the second disc; radial extent of the second disc; relative inclination of the discs; mass range and IMF slope 
of disc stars; binary fraction; simulated time; mass of the cusp SBHs; radial extent of SBH cusp; analytic cusp potential; post-Newtonian 
treatment for disc stars. All masses are i n sola r masses, all distances in parsec. Where present, both SBH and analytic cusps follow the 
density profile derived bv lSchodel et al. 1 1 |2007|) . For model 3d we assumed the multi-part power-law canonical IMF down to 0.01 Mq; to 
keep the total number of stars low, we had to replace the low-mass stars (< IMq) by IMq stars of equal total mass. The difference 
between the canonical and a = 2.3 models is that the latter are restricted to 1 — 120 Mq. 



0.3 (0.01 < m/M© < 0.08), on = 1.3 (0.08 < m/M© < 0.5), 
and «2 = 2.3 (0.5 < m/M©). The number of sta rs in the 
mass interval m to m + dm is £(m)dm, and the ISalpeteJ 
(| 19551 ) IMF has a = 2.35. Since the disc mass for an IMF 
as flat as a — 1.35 is not very sensitive to the lower stellar 
mass limit, we use a mass range of 1 — 120 Mq. 

Fo r most of our mode ls, we assume a disc age of 
5Myr (|Paumard et al.l [2006). The discs cannot be much 
older: According t o standard stellar evolution models (e.g. 
iHurlev et al.l |2000T ) , an age of 6 Myr implies a maximum 
zero-age main sequence (ZAMS) mass of 32 Mq for observed 
stars. The observed nu mber of disc stars wit h minimum 
ZAMS masses of 20 M (|Paumard et alj|2006h would thus 
imply a total initial disc mass of 16 750 Mq for the IMF with 
a = 1.35 (77 stars with 20 < m/M Q < 32), plus some bi- 
nary fraction, and the discs are belie ved to have a total mass 
<15 000M© |Navakshin et aLllioOtjh . However, we vary the 
disc masses and ages consistently in Sections 14.31 and 14.41 

To test for disruption of binaries, we mark a certain frac- 
tion /bin of stars in our models as binaries (varying between 
the different models, see Table [TJ, and double their mass 
(assuming equal-mass binaries). We track close encounters 
during our simulations to determine possible points of dis- 
ruption. Disruption of a binary of mass m,bm and separation 
a occurs if it passes another particle (SMBH, disc or cusp 
star) of mass M at a distance less than R tid ~ a yMjmun 
|Yu fc Tremaineil2003l ). 

We accounted for the influence of a background cusp as 
discussed in Section [2] either by using an analytic potential, 
or by adding 14 000 stellar black holes (SBHs) of 15 M Q into 
th e inner 0.22 pc, consis tent with the density profile derived 
bv lSchodel et alj (|2007h . 

Table [T] lists the varying parameters of the models 
used throughout this paper. On 3 GHz desktop comput- 
ers equipped with GRAPE6-Pro/8 or MicroGRAPE boards, 



calculation of these models took from five hours (top-heavy 
IMF, no SBH cusp) to a few 100 hours (all other models) 
per model. 



4 RESULTS 

4.1 Evolution of a single disc 

ICuadra et al.l (120081 ) have shown that a single cold disc 
in the absence of a perturbing potential cannot explain 
the observed large inclinations and eccentricities. The 
required perturbing pot ential may result fro m a mas- 
sive circum-nuclear disc l|Genzel et al.|[l994l ; cf. ISubr et"ai] 
2008), an intermed iate-mass black hole and /or star cluster 
(probably IRS 13; Maillard et al.l 12004 ), or a second d isc 
jpaumard et al.ll2006l ; cf. lLockmann fc Baumgardtj|2009h . 

To test whether a spherical cusp alone significantly al- 
ters the evolution of a disc, we integrated m odels of a single 
disc with the a = 1.35 IMF suggested by iPaumard et al] 



(|2006f ). restricted to 1 - 120 M©0 

Figure [1] shows the time evolution of eccentricity and 
inclination of a single disc in a s tellar cusp (la). In line with 
the results of lPerets et al.l (|200&f ) . we find a faster increase of 
eccentricities and inclinations than in the absence of a stel- 
lar cusp (lc). This still holds if we consider only the massive 
stars. However, less than 20 per cent of the massive stars 
(above 20 Mq) have final eccentricities e > 0.5, and a mere 
two per cent of the massive stars reach an inclination larger 
than 45° , while half of the disc stars observed have eccentric- 
ities e > 0.5, and a large number of stars have much higher 



3 Using a Salpeter/Kroupa IMF l lKroupall200lh with a = 2.3 for 
respective masses lead to similar results. For a discussion on the 
choice of underlying IMF see Section 14.31 
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Time / Myr Time / Myr 



Figure 1. Time evolution of orbital eccentricity (top panels) and inclination (bottom panels) of an initially cold disc. The left panels 
show the root mean square values for models in the presence (la; black) and absence (lc; green) of a stellar cusp. The corresponding 
dashed and dotted lines in the right panels show the fraction of stars on orbits with eccentricities larger than 0.3 and 0.5 (top), and 
inclinations larger than 20° and 45° relative to the fitted disc plane (bottom), respectively. While the thin lines refer t o the distribution 
of all stars, the thick lines represent the massive stars (> 20 Mq) only. In line with the results of lPerets et al.l d2008h . we find a faster 
increase of eccentricities and inclinations in the presence of a st ellar cusp. For comp arison, the blue points indicate the corresponding 
values for the young disc stars observed in the Galactic Centre llPaumard et al,ll2006l) , which have masses > 20 Mq. We assign Poisson 
error bars to the eccentricity values. Since no estimates are provided for the orbital planes, we use disc fitting of the velocity vectors for 
a very pessimistic lower limit of the inclination values. 



inclin ations to the observed clockwise disc ijPaumard et al.l 
2006). In fact, the largest contribution to the high values 
of observed eccentricities and inclinations stems from the 
counter-clockwise stars. 

We conclude that the stellar cusp does have a signifi- 
cant effect on the evolution of ecce ntricities and e speci ally 
orbital inclinations, but agree with ICuadra et al.l (|2008h in 
that it cannot account for the highly eccentric and counter- 
clockwise stars assuming an initial configuration of a single 
cold disc in the absence of any perturbing potential as dis- 
cussed above. 



4.2 Evolution of two interacting discs 

iNavakshin et al.1 (|2006l ) have shown how the mutual inter- 
action of two discs with large relative inclination leads to an 
increase of disc thickness. The evolution of orbital eccentric- 
ities and inclination (i.e. disc war ping) in the presence of a 
cusp h as been studied in detail by lLockmann fc Baumgardt] 
(2009). They find that while the evolution of orbital incli- 
nations is dominated by the discs' potential (precession and 
warping), the change in eccentricity and semi-major axis of 



the disc stars is driven by non-resonant as well as resonant 
relaxation with disc and cusp stars. 

Here we analyse the effect of the stellar cusp 
on the discs by computing the discs modelled by 
lLockmann fe Baumg ardt (2009) a) in a cusp of 15 Mq SBHs, 
b) in an analytic cusp, and c) in the absence of any cusp po- 
tential. 

Figure [2] shows the orientation of the orbital planes of 
the disc stars after an integration time of 5 Myr. Since pre- 
cession is dominated by the disc potentials, direction and 
amount of warping do not depend significantly on the cusp. 
On the other hand, the distribution of eccentricities does 
depend on the surrounding cusp: Figure [3] shows the fi- 
nal distribution of the disc stars' eccentricities as a func- 
tion of semi-major axis, for both single-disc and two-disc 
setup. While the overall eccentricity distribution is simi- 
lar in all cases, eccentricities e > 0.8 are only achieved in 
the presence of a stellar cusp, or in a two disc setup with- 
out any cusp potential (at the very inner edge). As we will 
discuss in Section 14.41 while high eccentricities reached by 
long-term Kozai interaction with the two discs are impeded 
by an analytic cusp, relaxation p rocesses in a stellar cusp 
(cf. lLockmann fc Baumga rdt 2009) can lead to eccentricities 
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Figure 2. Aitoff map of the orientation of orbital planes after an 
integration time of 5 Myr. The panels show the results of a model 
including stellar cusp, analytic cusp, and no cusp (top to bottom). 
Filled circles and squares denote the normalised angular momen- 
tum vectors of stars from the clockwise and counter-clockwise 
disc, respectively. The symbols are colour-coded by semi-major 
axis a; the open (red) symbols show the mean angular momentum 
vector of the initial discs. The direction and amount of warping 
does not depend significantly on the presence of a cusp. 



high enough for short-term resonances to be effective despite 
the spherical potential, thus explaining the high eccentrici- 
ties in three of our six cases (two discs in the absence of a 
cusp, or one or two discs interacting with cusp remnants). 



4.3 Constraining the disc mass and IMF 

As a direct re sult of the disc warping described in the pre- 
vious section, iNavakshin et all l|2006l ) have provided upper 
mass limits to the discs' m asses considerin g thei r appar- 
ent flatness as measured by IPaumard et alj (l2006t) after a 
few Myr of interaction. However. iLu et al.l l|2009h did not 
find convincing evidence for a counter-clockwise disc as flat 
as 19°, which can also be interpreted as an indication for 



a higher mass of the m a ssive disc in a two-disc model. 
lLockmann fc Baumgardtl l|200fjf ) have extended the argu- 
ment of warping as a mere measure of disc thickening by 
measuring the warping angle as a function of radial distance 
and found good agreement of a two-disc model with the ob- 
servations, showing how the warping angle can reach 50° 
already at 0.15 pc from the SMBH. 

Here we calculated different models with two discs of 
stars following a canonical IMF from 1 — 120 Mq. Using bi- 
nary fractions of 40 and 13 per cent implies total masses of 
27 000 Mq (3a) and 21 500 M (3b) , respectively, to have the 
same number of massive stars. Figure [4] shows the orienta- 
tion of orbital planes of model 3a after 5 Myr of integration, 
which shows the same trend as the results of the previous 
section; the same holds for model 3b. Assu ming a disc age 
of 4 M yr (the lower limit as estimated by IPaumard et al.l 
2006) reduces the initial disc mass by 20 per cent (3c), since 
a larger fraction of massive stars survives the shorter disc 
age, and thus a smaller initial number of massive stars is 
needed to explain the observed 77 stars. In another model 
(3d), we choose the initial disc m ass to be consistent with 
the canonical IMF |Kroupdl200lT ) down to 0.01 M , result- 
ing in 80 per cent more mass than with a cut-off at 1 Mq . 
To keep the total number of stars low, we had to replace the 
low-mass stars (< 1 Mq) by 1 Mq stars of equal total mass. 
Again, the overall spiral trend of disc warping as plotted in 
Fig. [4] is similar to the earlier results. 

Figure [S] shows the angle between a disc fitted to mas- 
sive stars at a certain radius and the outermost fitted 
disc. While the amount of warping of all clockwise (mas- 
sive) discs modelled with a canonical IMF is compatible 
with the observed values (or e ven below; cf. discussion in 
lLockmann fe Baumgardtl I2009T ) , the integration results of 
the counter-clockwise disc fit the observations well only out- 
side « 0.16 pc from the SMBH. Inside this distance, the 
counter-clockwise disc dissolves, and thus the warping angle 
is no longer a smooth function of central distance. However, 
it is essential to note that the observed stars are assigned 
to one of the two discs depending on whether they are mov- 
ing clockwise or counter-clockwise in the sky, since there is 
no other way to tell their disc of origin. This implies that 
none of the squares on the northern hemisphere of Fig. [4] 
would be considered for determining the warping angles 
of the counter-clockwise disc. Furthermore, more stars ob- 
served within 0.16 pc from Sgr A* move on projected radial 
orbits (and are thus completely dismissed) than on projected 
counter-clockwise orbits. If these (and even some of the 
clockwise-moving stars) stem from the more strongly warped 
counter-clockwise disc, they might significantly increase the 
amount of warping of the observed counter-clockwise disc. 
We waive a further analysis since assignment of stars to one 
of the two discs is rather uncertain in the inner 0.16 pc with 
the currently available observational data. In addition, the 
choice of the line of sight in our models is arbitrary. 

Considering this and the fact that there is as yet no 
model explaining all of the observed features of the young 
stars and including all perturbations like e.g. the CND, we 
conclude that these models are still consistent with the ob- 
servations. This seems to be in contrast to the results of 
INavakshin et al.1 l|2006l ) , however they used stellar discs only 
out to 0.2 pc from the SMBH, whil e more than one third of 
the stars listed by IPaumard et al.l l|2006l ) are further away 
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Figure 3. Eccentricity distribution as a function of semi-major axis after an integration time of 5Myr. The panels show the results of 
models of one disc (top) and two discs (bottom) including stellar cusp, analytic cusp, and no cusp (left to right). The stars are grouped 
by final mass above 15 Mq (black triangles), below 2Mg (open red squares), and in between (green circles). High eccentricities are only 
achieved in the presence of a stellar cusp, or by interaction between the two discs in the absence of a spherical potential. 
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Figure 4. Orientation of orbital planes of model 3a following 
a Kroupa IMF above 1 Mq after an integration time of 5 Myr. 
Despite significantly higher disc masses, the overall trend of disc 
warping is comparable to the results displayed in Fig. [2] 



from Sgr A*. In addition, they did not distinguish between 
stars in the outer parts of the discs which almost retain their 
orbital planes, and stars close to the centre which strongly 
precess. 

Our results show that a canonical IMF over the range 
0.01 — 120 Mq cannot be exclu ded dynamically for th e young 
discs in the Galactic Centre. IPaumard et alJ l|2006l ) find a 
flat IMF from analysing the if -band luminosity function of 
the disc stars; however, they assume a lower mass limit of 



20 Mq for the observed stars, and their estimate of the disc 
age of ~ 6 Myr limits the initial masses of the now observed 
stars to ~ 32 Mq, so this slope ref ers to a fairly small mass 
interval only. iManess et al.l ([2007) find a best fit of obser- 
vations in the central parsec of our Galaxy with a model of 
conti nuous star formation w i th a top-heavy IMF. Further- 
more, [N^^kshi^^Sunyae^ (H005|) argue that the X-ray lu- 
minosity of the Sgr A* field is one order of magnitude too 
low to account for the number of young < 3 Mq stars ex- 
pected from a canonical IMF, considering the large number 
of O-stars observed in the discs. This can also be explained 
by a higher low-mass c utoff of the IMF at of the order of 
1 Mq . [Alexand er et al.l (|2007T ) find that eccentricities as high 
as 0.2 — 0.3 can only be explained by a top-heavy IMF when 
assuming a circular disc origin, however they neglect the 
influence of the cusp and the s econd disc as elaborated by 
lLockmann fc Baumgardt! ||2009l ); both suggest an eccentric 
disc origin for the counter-clockwise system. Simulated star 
formation of gaseous discs suggests that stars form with a 
top-heavy IMF only if the cooling times ar e .just as large as 
allow ed to still induce disc fragmentation (|Navakshin et all 
120071 and references therein) . Unfortunately, theoretical IMF 
predictions have failed in the past to c orrectly descr ibe the 
observations near the Galactic Centre (Kroupa 2008). 

Another argument for a top-heavy IMF arises from the 
relation between disc and cusp mass. Assuming that we do 
not live in a special time, the formation of stellar discs in 
the Galactic Centre must be a frequent event. To explain the 
observed number of > 20 Mq disc stars at an assumed age of 
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Figure 5. Disc warping angle as a function of central distance. 
We plot the angle between the plane fitted to any disc's subset 
of stars and the outermost fitted plane for the respective disc. 
The left panels show the results for the clockwise disc from ob- 
servations (grey error bars) as well as for models 3a through 3d 
(blue lines, top to bottom); the right panels show the respective 
results for the counter-clockwise systems. Horizontal and verti- 
cal error-bars indicate standard deviation of central distance and 
f cr-thickness of the fitted planes, respectively. The corresponding 
error bars for the simulations are generally smalle r and thus omit- 
ted for clarity (cf . lLockmann fc B aumgardt 2009j) . In none of our 
models does the amount of warping of the clockwise disc exceed 
the observed value. The modelled counter-clockwise discs dissolve 
towards the inner part, where in most cases warping is no longer 
a smooth function of central distance. 



5 Myr with the a — 1.35 IMF suggested bv lPaumard et all 
([20063). a total initial disc mass of ~ 15 000 Mq is required, 
of which a mass of ~ 3000 Mq ends up in stellar remnants. 
Using the canonical IMF requires a total initial disc mass of 
~ 45 000 M , of which ~ 30 000 M in remnants and low- 
mass stars are retained over a Hubble time. Assuming that 
most stars survive in the central region (or are replaced by 
inspiralling massive remnants), the observed cusp mas s of 
~ 750 000 M in the central 0.5 pc (|Schodel et alj|2007f ) al- 
lows for the formation of discs like those observed every 
50 Myr for a top-heavy IMF, and only every 500 Myr for a 
canonical IMF. The fact that we observe very young stel- 
lar discs thus seems much more unlikely if star formation in 
the Galactic Centre follows the canonical IMF. However, it 



may also be an indication of recently enhanced star forma- 
tion processes. For example, if the Galactic Bar is young, as 
an increa sed fraction of ba rred galaxies for lower redshifts 
suggests |Sheth et al.ll2008f ) . bar-induced gas inflow may ex- 
plain such an enhancement |Sellwood fc Wilkinsonlll993l and 
references therein). 

We conclude that while theoretical star formation ar- 
guments and some observations indicate a top-heavy IMF 
in the centre of the Milky Way, our self-consistent models 
of stellar dynamics allow the existence of a canonical IMF 
as well. The observational evidence clearly needs to be im- 
proved by star counts, while the theory of star formation 
is currently not able to make conclusive statements o n the 
IMF of stars in the Galactic Centre (see lKroupall2008l for a 
discussion of recent failures of IMF theory) . 



4.4 Creation of S-stars in a cusp environment 

In lLockmann et all (|2008l ) we have shown that a signifi- 
cant number of disc stars eventually achieve highly eccen- 
tric orbits due to Kozai interaction, thus providing a nat- 
ural mechanism to create the observed S-stars and hyper- 
veloc ity stars by d i srupt ion of close-passage binaries. How- 
ever, Ivanov et alJ (|2005h remarked that the effects of peri- 
centre shift due to relativistic precession and a spherical 
cusp can destroy the Kozai mechanism. While the relativis- 
tic precession period decreases strongly with increasing or- 
bital eccentricity, 



. GH = 2 X 10°P orb 



0.1 pc 



(1-e 2 ) 



(2) 



(where P rt>, a and e are the orbital period, semi-major axis, 
and eccentricity, respectively), they showed that cusp pre- 
cession at a given distance is most effective for near-circular 
orbits: 



73Por 



0.1 pc 



1 



(3) 



(both equations assuming Msmbh = 3.5 x 1O 6 M and the 
cusp profile inside 0.22 pc bv lSchodel et al.ll2007l ). From this 
we conclude that the very high eccentricities in Kozai reso- 
nance are prevented by relativistic precession at high eccen- 
tricities and by cusp precession at low eccentricities. 

Figure [6] shows the eccentricity evolution during Kozai 
resonance for integrations of a three-body problem with and 
without a cusp, and with and without PN treatment for rel- 
ativistic effects. It is clearly seen that the PN terms prevent 
higher eccentricities when a certain eccentricity threshold 
is reached (in this case e ~ 0.9), while their effect on the 
rest of the Kozai cycle is not dramatic Q This can be un- 
derstood by general relativity being most effective for high 
eccentricities (and thus small pericentre distances). On the 
other hand, the changes on both maximum eccentricity and 
Kozai period exerted by a stellar cusp are quite low, un- 
til a certain mass threshold of the cusp is reached, which 
prevents not only the highest eccentricities but essentially 



4 The minor changes in oscillation period caused both by cusp 
and PN terms reflect the small perturbations to the Kozai cycle 
due to precession, where resonance is slightly out of phase with 
the unperturbed case. 
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Figure 6. Eccentricity evolution of a test star undergoing Kozai 
resonance driven by a fictious 15 000Mq particle representing 
the disc potential. Both particles have initially circular orbits 
about the 3.5 X 10 6 Mq SMBH with semi-major axes of 0.04 pc 
and 0.16 pc, respectively, resulting in the star's orbital period of 
400 yr. We calculate models without cusp (solid black line), one 
per cent (dashed line), and five per cent (dotted line) of the ex- 
tended cusp mass observed in the Galactic Centre (modelled as a 
smooth potential). Each curve is accompanied by a corresponding 
thick green curve which represents a respective integration includ- 
ing post-Newtonian terms up to 2.5 PN to account for the effects 
of general relativity. While relativistic effects damp the Kozai ef- 
fect at high eccentricities, the presence of a stellar cusp leaves 
it almost unchanged until a mass threshold of a few per cent of 
the observed value, where eccentricity growth is stopped at very 
small values so that the orbit remains almost circular (straight 
dotted lines). 



any significant eccentricity increase and thus the complete 
Kozai cycle. Even for the innermost disc stars (a ~ 0.03 pc), 
we find this thresho ld to be at on l y a fe w per cent of the 
cusp mass given by ISchodel et al. | {20071 ): While the 1 per 
cent curve reaches e = 0.99, the 5 per cent curve does not 
exceed e = 0.01. 

As a result of this analysis, we cannot expect any disc 
star in the Galactic Centre to undergo a full Kozai cy- 
cle, which was also shown by IChand (|2009l ). However, if a 
star achieves high e ccentricity by (non-)resonant re laxation 
or close encounters (|Lockmann fc Baumgardtll2009l ). it may 
find itself in a point of Kozai resonance where the time until 
maximum eccentricity is significantly shorter than the full 
Kozai period, and also the cusp precession period is longer 
than for the circular case. Furthermore, in our system of two 
disc potentials and a large reservoir of perturbing disc and 
cusp stars, other forms of resonance may occur (short-term 
Kozai/resonant relaxation with changing sources of pertur- 
bation) . 

Going beyond lLockmann et all {2008), we computed 
different models of two-disc systems in the presence of a 
stellar cusp. As a result of cusp precession, we find only 
1 — 2 binaries in each model achieving eccentricities high 
enough for binary disruption to occur (4a, 4b). When using 
a canonical IMF down to 1M©, 7 S-stars are created (4c). 
Besides the smaller disc extent used for model 4c (resulting 
in faster relaxation and shorter Kozai periods), this is due 
to the number of B-stars being four times higher than in the 
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Figure 7. Evolution of orbital eccentricity (top) and x compo- 
nent of the normalised angular momentum vector, jx/j (bottom), 
of three selected stars of model 4b. The cyclic nature of eccen- 
tricity evolution as well as the rapid change of the stars' orbital 
planes as indicated by jx/j are revealing of the Kozai mechanism. 



previous cases (and also the total mass increase), while the 
number of O-stars is still consistent with the observed value. 
Due to the increase of total disc mass, the number of S-stars 
created may increase even further when using a canonical 
IMF down to 0.01 M Q . 

For those stars achieving high eccentricities we can see 
the behaviour as predicted above: The eccentricity changes 
in a random fashion by relaxation processes until a certain 
value of high eccentricity is achieved. At this point, the star 
(or binary) goes into Kozai resonance for a time scale long 
enough to achieve extreme eccentricity, but too short to be 
prevented by precession effects. 

Figure [7] shows the evolution of eccentricity as well as 
one component of the normalised angular momentum vec- 
tor, j x /j, of selected stars of model 4b achieving high eccen- 
tricity. While the eccentricity evolution seems at first sight 
compatible with the secular instability in flat eccentric discs 
as discovered bv lMadigan et all {2009), the rapid change of 
the stars' orbital planes as indicated by j x /j is revealing 
of Kozai resonance. The former secular instability can only 
occur in a flat (single, non- warping and non-thickening), ec- 
centric disc in a cusp with a density profile p <x r~ 15 . 

Interestingly, we find the orbital planes eventually flip- 
ping at very high eccentricities (e.g. solid curve at maxi- 
mum eccentricity, t m 5.1 Myr: j x /j suddenly changes from 
~ 1 t o ps — 1). This was first discovered by iTouma et al.l 
(2009) who applied their softened version of Gauss' method 
to the two-disc problem. However, in our case no softening 
is involved. At such high eccentricities, no large torque is 
needed to change the small angular momentum to the op- 
posite direction (note that this does not occur in an isolated 
three-body system undergoing Kozai resonance). 

In another calculation (4d; using the top-heavy IMF 
with a = 1.35), we used the PN terms for all disc stars 
from the beginning, and a spherical cusp composed of 40 Mq 
SBHs, which resulted in the creation of 8 B-type and 1 
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O-type S-star. Almost all of the originating binaries had 
close encounters with one of the SBHs (although not close 
enough to be disrupted), however none of these encounters 
lead to a sudden change in orbital eccentricity or semi-major 
axis which could explain the high eccentricity it achieved 
much later in the integration. We ascribe the higher num- 
ber of S-stars created to the fact that the effectiveness 
of resonant relaxation is proportional to the SBH mass 
ijRauch fc Tremainell 1996T ). Apparently, pericentre shift due 
to relativistic effects does not prevent such high eccentrici- 
ties. 

We conclude that if the stellar cusp in the Galactic Cen- 
tre is mass-domin ated by massive SBHs (as pred icted by 
iFreitag et al ] |2006l and I Alexander fc Hopmar]|2009r i , at least 
some of the S-stars observed could have originated from the 
young stellar discs. In a setup of an initially e ccentric disc (as 
suggested bv lLockmann fc Baumgardt) 2009 for the counter- 
clockwise disc), we can expect a somewhat higher probabil- 
ity of high-eccentricity orbits required for the formation of 
S-stars. Since the age of the S-stars is less well-determined 
than that of the stellar discs, some of them may stem from 
an earlier generation of discs, which would anyway be needed 
to explain the longer travel times of the hyper-veloc ity stars 
observed in the Milky Way halo l|Brown et al.|[2007l and ref- 
erences ther ein), assuming the sam e disruption scenario as 
their origin l|Lockmann et alT 2008). Furthermore, scattering 
of young binaries by massive pe rturbers may cont ribute to 
the observed number of S-stars l|Perets et al.ll2007T ). 



5 CONCLUSIONS 

In this paper, we have investigated the influence of the stellar 
cusp on the evolution of young stellar discs in the Galactic 
Centre. Our results show that the impact of stars and stel- 
lar remnants on the disc stars' orbits must not be neglected 
in theoretical models and numerical simulations: The stel- 
lar cusp enhances evolution of orbital eccentricities and disc 
thickness by (non-)re sonant relaxation of angu l ar mo men- 
tum as described by lLockmann fc Baumgardtl (|2009l ). We 
find that this is especially true for high eccentricities, which 
are prevented by a smooth cusp potential resulting in orbital 
precession, but are achieved in a grainy cusp as we expect to 
find it in the inner parsec of our Galaxy. A combination of 
two-body relaxation processes resulting in a random increase 
of eccentricity and short-term Kozai resonances leads to ec- 
centricities high enough for binary disruption, leaving be- 
hind an S-star on a tight orbit and ejecting a hyper- velocity 
star. This process is not as efficient as it would be in the 
absence of a stellar cusp, but a few generations of stellar 
discs could account for the observed S-stars. 

We conclude that the creation of S-stars is a natural 
consequence of the interaction of two stellar discs in a grainy 
cusp around an SMBH, where stellar orbits are exposed 
to the effects of general relativity. Scattering of stars onto 
highly eccentric orbits is enhanced by a cusp composed of 
more massive particles, such as SBH multiples, (remnants 
of) star clusters or IMBHs. Thus, more realistic modelling 
of stellar dynamics in the Galactic Centre requires a better 
understanding of the composition of the stellar cusp. 

Even more so, the composition of the stellar discs re- 
mains unclear. The mere existence of the S-stars may be 



regarded as evidence for a steeper (i.e. more standard-like) 
IMF, as this would significantly increase the number of B- 
stars (and thus stars eventually reaching tight orbits) given 
the observed number of O-stars. A detailed survey of B-type 
stars in the central parsec of the Milky Way will not only 
help us to better understand the composition of the young 
stellar discs and thus the formation of the S-stars, but also 
increase the significance of the counter-clockwise disc as well 
as the observed disc warping. However, if the formation of 
stellar discs in the Galactic Centre is a common scenario, 
most of the B-type stars (including the known S-stars) may 
stem from previous generations of discs, and their orbits are 
expected to be randomised by the influence of the cusp as 
well as the potential of current and previous discs. 

As a by-product of our models, we can determine the 
merger rates of stellar black holes with the SMBH. These 
mergers can be caused by an SBH plunging into the SMBH, 
or by gradual inspiral due to gravitational wave emission 
(so-called extreme mass ratio inspirals, EMRIs), yielding a 
promising scenario for the detection of gravitational wave 
signal s by the planned Laser In terferometer Space Antenna 
(e.g. lAmaro-Seoane et al.M2007f ). Since in our models we fo- 
cus on the evolution of the stellar discs, we do not differenti- 
ate between EMRIs and plunges of SBHs, but only identify 
the merger of an SBH with the SMBH as the point where 
their distance is smaller than three Schwarzschild radii. Our 
models including a cusp of 14 000 SBHs of mass 15 Mq 
showed merger rates from 20 — 80Myr _1 , with a mean of 
58Myr _1 . Note that these rates may be higher when using 
SBHs beyond 0.22 pc or of smaller masses (and thus a higher 
number). On the other hand, they need to be corrected for 
cusp depletion and loss-cone-refilling. 
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